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Abstract 



We present a calculation of the QCD correction factors r]2B and ri2K up to Next-to-Leading Order 
\l I within the MSSM. We took into account the region of low tan/? for the Higgs- and chargino sector 

I while neglecting the effect of gluinos and neutralinos. 
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^ ■ 1 Introduction 

Q^', During the last decade, Supersymmetry (SUSY) has been one of the most popular candidates 
for physics beyond the Standard Model (SM), see, e.g., the reviews 10-0]. Among all other ideas 
SUSY seemingly has the greatest potential to allow for the construction of a Grand Unified 
Theory. Even more, a low-energy version of SUSY would cancel the quadratic divergencies 
emerging within the Higgs-sector and cure the most nagging shortcomings of the Standard 
5^ I Model 0-0. However, so far there is no experimental evidence for any of the new particles 
predicted by the various supersymmetric models up to scales of the order of the W- and Z- 
mass. So before the advent of the next generation of new colliders one probably has to rely on 
indirect probes to test the idea of a low-energy SUSY. 

In this respect the most promising candidates are processes involving Flavour Ghanging Neutral 
Currents (FGNG), especially rare decays like 6 — > 37 and oscillations in the system of neutral 
mesons, namely BB— and i^i^— mixing [|TU|-[PU|. In this paper we will focus on the latter class 



of processes and consider contributions to the mixing phenomena mediated by the chargino- 
squark and the charged Higgs-up-type quark sector of the Minimal Supersymmetric Standard 
Model (MSSM) taking into account gluonic Next-to-Leading Order QGD-corrections (NLO). 
We present results for the mass-splitting AM^ within the corresponding neutral 5-system and 
for the parameter describing GF-violation in the i^-system. 

Within the framework of the SM the corresponding AF = 2 processes are governed by box- 
diagrams with internal up-type quarks and W^-bosons and are suppressed by either small quark 
masses {u and c) when compared to the l^-mass or by small mixings between the third and 



the first two generations. The corresponding Leading Order diagrams of the SM are depicted 
in Fig. ^ The result for the SM up to NLO QCD corrections can be found in [^-[^ . 

Within the MSSM - the supersymmetric extension of the SM with minimal field content, the 
same gauge group and soft symmetry breaking - there are different sources of flavour-changing 
processes. Nevertheless, all of them enter the stage of AF = 2 transitions via box-diagrams. 
First of all there are loops with internal up-type quarks and charged Higgs-Bosons (see Fig.|^). 
Taken alone, this corresponds to the Two Higgs Doublet Model (THDM), another popular 
extension of the SM. Results up to Next-to-Leading Order have been presented in pOl. The 



second source are loops involving charginos and up-type squarks (see Fig. y). Their effect 
is investigated in this article up to NLO of gluonic QCD corrections. The last source of the 
AF = 2 processes within the MSSM are off-diagonal elements within the squark-mass matrices, 
entering via boxes containing down-type squarks and gluinos or neutralinos (see again Fig.|^). 
In the context of this paper we do not take into account their effect. Additionally we do not 
consider supersymmetric Q CD-corrections mediated by gluinos although their effect might be 
sizable even in the limit of heavy gluinos [^. The incorporation of these contributions would 
require additional extensive calculations which is beyond the scope of the current paper. In 
Leading Order the effect of the various MSSM-contributions on the mixing phenomena has 
been examined in 



There it is shown, that i^'i^'-mixing strongly limits the size of some 
off-diagonal elements within the down-type squark matrix. Additionally, within the context of 

In this specific version 



the mSUGRA model ||3^-||3^, LO-results have been presented in |^0| 
of a supersymmetric extension of the SM, the soft supersymmetry-breaking terms are governed 
by just 5 parameters, reducing the vast number of parameters occuring within the MSSM 
considerably. The authors of claim, that the most sizable contributions to AMb and ex 
of the mSUGRA model in addition to the SM stem from the chargino-squark and the charged 
Higgs-up-type quark boxes when taking into account constraints of the parameter space as 
given by various other processes. These findings partly justify our focus on the charged Higgs- 
and the chargino-sector as primary sources of FCNC-processes. 

It should be stressed here, that a priori the ratio of the vacuum expectation values of the two 
Higgs doublets, tan/3, is a free parameter. Overall fits to available experimental data constrain 
this parameter to two regions, tan/5 ~ 30 — 50 and tan/5 ^1 — 2, and it seems that the latter 
one is the more favourable one pT| , Here we will concentrate on this region of low tan/?, 

which yields the major contributions to the mixing phenomena as stated by |^ . 



The paper is organized as follows. After we redisplay the necessary parts of the MSSM-Feynman 
rules we shortly review the basic formalism in LO. In section 3 we discuss some features of the 
NLO calculations with special focus on the matching-procedure. In section 4 we provide the 
results in NLO for the different observables and scan to some extent the parameter space of 
the MSSM. We close with some concluding remarks. 



2 Notation and basic formalism 
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2.1 Notation and Feynman— Rules 



Throughout this paper we use the notation of Here it should be sufficient to hst the 

relevant Feynman-rules and define the quantities involved. For the gauge bosons we employ 
the Feynman-t'Hooft gauge ^ = 1, the propagators can be found in the Appendix. The relevant 
vertices are depicted in |l|. There, the capital letters J, J, K denote the quark generations, the 
letters z,j label the squark-, Higgs- or chargino-fields and the a, 5, c are SU(3) indices in the 
appropriate representation. The matrices and 2^[7,d diagonalize the mass matrices X of the 
charginos and the up- and down-type squarks respectively, 

{Z-fX^Z^ = diag(mx, J , zIjjXu,dZu,d = diag(m„-^^^-J, (1) 

Throughout this paper we will assume a specific form of the matrix Z^ responsible for the 
mixing of the different squarks, namely all off-diagonal elements zero besides the ones associated 
with the mixing of the stops, m^'®. Additionally, we assume the up-type squarks of the first 
two generations degenerate in mass. In other words, we assume, that the squark mass-matrix 
preserves fiavour and mixes only the stops of different chiralities. This special form of the 
Z-matrices used in this paper is given in Eq. (|39D. It should be stressed here, that the last 
two assumptions concerning the up-type squark sector are mainly for the sake of a compact 
presentation but still somewhat unfounded. However, a generalization of our results to a more 
complicated structure of the up-type squark mass-matrix is quite straight-forward. It should 
be noted, that the four-squark vertex containing two up-type and two down-type squarks is 
not of relevance when considering 0(as)-corrections. 

C^"^ denotes the CKM-element between generations I and J and tan/5 = V2/V1 is the ratio of the 
vacuum expectation values of the corresponding Higgs-doublets. We have explicitly written 
down the Yukawa-type couplings in terms of the quark masses involved. Simple inspection 
shows, that for the region of low tan/5 we can neglect the couplings proportional to the down- 
type masses with respect to the top-mass. 
We use the short-hand notation 

Pl,r = , Y'^ = YPl (2) 

for the left- and right-handed projectors and their combinations with gamma-matrices. 



2.2 Leading Order 

We want to summarize now briefly the ingredients used to obtain Leading Order results. In case 
of the BB-mixing within the Standard Model this is quite transparent. Neglecting all masses 
besides the W- and the top-mass one can use the GIM-mechanism to account for the two lighter 
quark-types running in the boxes and calculate the remaining box-diagrams. This is simplified 
even more by using the Feynman-t'Hooft gauge for the W, leaving us with the physical W- 
boson and the would-be Goldstone-boson $, both with identical masses. Results for boxes 
involving additional Higgs-bosons and top-quarks can be obtained in the same manner. Due 
to the fact that the supersymmetric partners of the SM-particles differ in Spin by 1/2 one has 
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to perform a Fierz-transformation for diagrams with squark- and chargino-lines to write down 
the emerging operator in terms of colour-singlet currents. As stated before, we assume, that 
the squarks of the first two generations [ui ,i = 1, 2, 4, 5) are degenerate in mass and do not mix 
considerably. This allows us to use the GIM-mechanism within the squark sector to tighten 
our final results. Generalization to non-degenerate squarks is straightforward. 
The procedure of integrating out the heavy degrees of freedom leaves us with one operator 
and its Wilson-coefficient at the matching scale /xq = 0{M^) and yields the starting condition 
for the renormalization group equations. The connection between the starting scale and the 
bosonic scale fi = O(m^) is given in LO by the diagrams of Fig.^ and a corresponding factor rj2. 
Putting everything together, the effective Hamiltonian for AB = 2-transitions at the bosonic 
scale is given by 

^eff = 2 ^ m'w {VtdV,lf V2 S OlM , (3) 

where 

Oll = {da,Lk) {d,r%) ■ (4) 

We have written out the colour indices for the quark fields explicitly. The coefficient S" is a sum 
of the well-known Inami-Lim functions for the different internal particles, 

S = S{xw,xh) + S{{xi,ya}) 
S{xw,xh) = Sww{xw) +'^Swh{xw,xh) + Shh{xh) ■ 

2 6 _ _ 

S{{xi,ya}) = J2 Kij,abS{xi,Xj,ya,yb) (5) 

i,j=l a, 6=1 

Throughout this paper we use the following abbreviations 

^^^W,H ^"W ^^^W ^^^W 

for the ratios of masses entering the Inami-Lim functions. 

The K in Eq. (|^) account for the squark-quark-chargino couplings as given in the Feynman- 
rules above. They read 

K^,ab= n n U'Z^, ^ J" . Z\j^'''Z^\ . (7) 

c={a,b} k={ij} V V2Mw smp J 

For the i^-system the situation is slightly different. Because of the much lower mass of the 
fT-meson, the charm-mass can not be neglected any more, and one has to account for this fact 
by integrating out the heavy degrees of freedom in consecutive steps So we end with a 
modified effective Hamiltonian for AS* = 2-transitions at the bosonic scale, namely 



1 G 



2 

F _2 



4 TT 



\Vc,v:,)\ Six,) + {Vt,V*)% S 

+iVcdV:,)iVt,V*)v, Six,, xw)] OlM , (8) 
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where the contribution of the heavy particles alone, S is given by Eq. (^) and the contributions 
involving a c-quark are contained in S{xc) and S{xc, xw) with corresponding QCD-correction 
factors ?7i,3. Of course the operator Oll changes accordingly, the 6-quarks are replaced by 
strange quark fields. 
The factors rj2{B,K} in LO are given by 



ri2B,LO = OisiMw 



V2K,LO = asimcy '^''-^ M — — r — - — - , (9) 

\as[mc) ) \ asimb) J 

where we have explicitly chosen the matching scale /xq = Mw- 'J^^^ and Z?^'^-' are the coefficents 
of the anomalous dimension of the operator Oll and the QCD /3-function in LO, the latter 
one is labelled by the number of active flavours, nf. They are given by 

,10) 

Nc = 3 denotes the number of colours. We will see in the following, that we are able to plug the 
NLO QCD-corrections within the MSSM into the corresponding factor ri2, see Eqs. (0,1). The 
other QCD-factors rji and 773 remain unaltered. Because they are even in LO quite complicated 
functions of a^, we just quote the numerical results for the NLO expressions we will use later, 

rii ^ L38 , ri3 ^ 0.47 (11) 

as given in 



The last task left is to extract the physically observable quantities from the effective Hamilto- 
nians. In the case of the mass difference AMb in the i?-system and ex one has to sandwich 
the Hamiltonian in between two mesonic states. This amounts to the use of the relation 



{B\OLL{f^)\B) = ^BBifi)fBmB, (12) 



where the so-called bag-parameter is given scale-independently by 



Bb = 5B(/.)a.(/.)^'"/(<) 



1-^Z 
An 



(13) 



Similar expressions hold for the i^T-system as well. In LO the last term of equation (|r^) can 
be omitted, the term appearing there will be defined in the next section. 
Unfortunately the evaluation of and BbJb respectively by lattice calculations yields large 
uncertainties 



Bk = 0.75 ± 0.15 , ^BbJb = 0.2 ± 0.04 GeV , (14) 

and it is fair enough to state that this fact spoils the proper determination of the CKM-element 
Vtd by the processes considered here. 
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3 Explicit QCD— corrections, matching and running 



3.1 Explicit QCD— corrections 



We will now proceed by presenting the explicit perturbative QCD corrections up to 0{as)- 
They are obtained by evaluating the diagrams in Figs. ^ and ^ For the diagrams containing 
quarks and bosons we performed the calculation in an arbitrary covariant ^-gauge for the gluon 
and - as stated before - the Feynman-'tHooft gauge for the VT-boson. For diagrams involving 
supersymmetric particles we chose explicitly = 1 for the gluon. Again, we perform a Fierz- 
transformation for the latter diagrams to write down the operators in terms of colour-singlett 
currents. 

As one easily notices, diagrams (b,c,f,g,h) have the octet-structure Ta®T°- whereas the diagrams 
(a,d,e,j) have the singlet-structure 1® 1. The double penguin diagrams (h) and the diagram (i) 
involving the four-squark coupling in a similar fashion do not contribute for vanishing external 
momenta. 

Furthermore, we have to face ultraviolet as well as infrared divergencies. The first ones stem 
from the diagrams (d,e,j) . In this case we employ dimensional regularization within the so- 
called Naive Dimensional Regularization-scheme (NDR-scheme) and we renormalize 



using the M5'-scheme ^\ Note, that for the supersymmetric particles we have to encounter 
additionally tadpole-like contribution (i) caused by the four-squark vertex. The infrared diver- 
gencies stem from the diagrams (a,b,c) . To handle them, we keep the masses of the external 
quarks whenever necessary. We will see, that this particular choice does not affect the final 
result for the Wilson-coefficient after performing the matching. Of course we could have used 
dimensional regularization as well, but the method chosen allows to compare the calculation 
presented here step by step with the ones presented earlier PU[ . 

In the following we will only consider the contribution to 772, the parts of (H) involving the 
charm quark will not be displayed in this section. This reduces to (^, and we will perform 
the necessary steps to gain rj2B at NLO in this section before we just generalize on r]2K at the 
end. 

So the 0(as)-corrections to the Hamiltonian of (13) show the structure 



AH, 



eff 



47r2 

with U given by 



(15) 



(16) 



where a sum over = 1, 2, 3, LL is performed. The operators Ok read 



Oll 
Oi 
O2 
O2 



{L^R}, 
+ {L^R}. 



(17) 
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The operators Oi,2,3 stem from the infrared divergent diagrams b,c and a, respectively. All of 
them are written in a Fierz-invariant fashion, allowing us to perform the Fierz-transformation 
for the supersymmetric contributions. 

The colour factor Cp along with Ca to be used later is given by 
A^^ - 1 ~ - 1 

^^^^aT'^^^^^- ^^^^ 

In analogy to Eq. the functions can be decomposed as 



(19) 



for the parts involving quarks and bosons and the parts with squarks and charginos, respectively. 
The functions x^j^^i^w, xh) have been calculated already 0, here we just restate the result 



Xll{Xw,Xh) 



X^l{xw,x,) 

Xf\xw.xn) 
Xfix^.xn) 



+ 



2i{l + giR + ln(a;^J) + 6 ln(a;^J ^ 



Xjd 
dxi 



S{xw, xh) , 



L^''™^'\xw,xH) 

+ [2^ (1 + qir) + (3 + In(xbXd)] S{xw, xh) 
-{?> + i)S{xw,XH) , 



ml - ml 



Xb 



Xd 



In — S{xw,xh) ■ 



All other x equal zero. We have introduced here some other abbreviations, 



X 



MO 



9lR 



Xd\ia.{xd) - Xfcln(xb) 



(20) 



:2ii 



where fio is the scale of integrating out the heavy degrees of freedom. Note, that all masses 
entering the functions are the masses taken at this scale, for example, mt is the top-mass 
renormalized within the MS'-scheme at the scale /zq- 

For the squarks and charginos running in the boxes the functions x differ only slightly from 
the X just presented. Namely, for the Xi,2,3 one just has to replace S{xw,xh) by S{{xi,ya})- 
The xll read for ^ = 1 



XLLii^t^ya}) 



2 6 _ ( 

^iho,b \ L^^\x^j, ya,b) + [2(1 + giR + In(x^J) 

i,j=l a,b=l \ 

+ 21n(x,JE E (^ + 

k=i,j c=a,b \ * 



'Sxid ^ 2ycd' 
dye , 



S{Xij, ya,b) 



XLL{{Xhya}) 



Yl Y Kij^ab{L^^\xij,ya,b) 
i,j=l a,b=l 

+ [2{l + giR)+Aln{xbXd)]} S{x 

ya,\ 



(22) 
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To perform the calculations we used Mathematica [^] and the Mathematica-based package 
FeynArts |50| . 



Note that the results so far show a non-trivial dependence on the gauge parameter and the 
IR-masses. This dependence will vanish after we matched the full and effective sides of the 
theory. The functions L constitute an integral part of the Wilson-coefficient emerging after 
this procedure, and are gauge and IR-independent. They are listed in the Appendix. 



3.2 Matching and running 

To match the full and the effective side of the theory we have to evaluate the matrix element 
of the physical operator On up to order 0{as) using the same regularization, renormalization 
and gauge prescriptions employed before. Yet, the one-loop amplitude of Oll as given by Fig.^ 
results in the same operators with the same coefficients for the obviously unphysical operators 

^1,2,3, 

Hoop tree 

E {Cpxi]M + CAxtM) 1 ® 1 (^fc)t- , (23) 

where again the sum over k = LL, 1, 2, 3 is understood. We have used the identity 

Ta^T" = Ca1®1 (24) 

to allow for rearranging the colours under Fierz-transformations. Now the Xa(/Uo) read 

Xa,lM = -3 + 2^{l + \n{xJ+ gin), 
Xa?ll(/^o) = -5 + 2^(1 + gin) + {3 + ^)\n{xbXd)-Qln{Xf,g), 
x2!i(/io) = -(3 + 0, 

XA.sif^o) = ^Xa,2(/^o) = ^-^^J^^^ln f— ) , (25) 
2 ' 2 - \XdJ 

and equation (^) enables us to write the matrix element of the effective Hamiltonian up to 
order 0{as) as 

(i/eff) = {Hi^ + AHii) = ^M^{V,dV*yCLM{OLM)iioor>. (26) 
We are now able to expand the Wilson-coefficient Cll to 0{as)^ 

CLLif^o) = S+^D (27) 



with S given by equation (|]). Comparing coefficients yields the result for D, 



D 



D{xw, Xh, X^o) + D{{Xi, Ua}, x^J , 



3 + 61n(x^J Y 



Xjd 
dxi 



S{xw,xh) 



D{{xi,ya},x^ 



+Ca{l^''^''''^'\xw,xH) + [5 + 61n(x^J] S{xw,xh)} 

2 6 _ ( 

Cf J2 Y Kij,ablL'^^Hxij,ya,b) 
i,j=l a,b=l (. 

+Ga J2 Y1 ^iJM \ L^^\xi,j^ya,b) + [5 + 61n(x^J] S{xij,ya,b) \ ■ 

i,j=l a,b=l I J 



S i-^itjy ya,b) 



(28) 

Hence, at this stage we have expressions for the Wilson-coefficients at the matching scale /iq, 
the last step to be performed is the evolution down to the mesonic scales. 
This evolution is achieved with help of the renormalization group. The solution of the renor- 
malization group equation for the Wilson-coefficient Cn at NLO reads 



Cll{^) = V2,LO 



An 



w 



(29) 



with the matching scale chosen explicitly as /xq = My/. The term is given along with other 
ingredients needed as 



^(1) 

Irif 



iVc-1 



-21 



2Cfnf 
57 19 



Nc + -Uf 

Nr 3 3 \ 
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(0) 



2Mr' '"' 



(30) 



The leading order expressions have been given in Eqs. (0, 110). Defining 



O 



LL = a. 



(31) 



and setting nj = 5 for the 5-system we can cast the effective Hamiltonian for I\B = 2- 
transitions at the mesonic scale in the following form 



H, 



eflf 



47r2 



M^{VtdV,l)%SOLL 



(32) 
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where the factor 772 describing QCD-corrections up to Next-to-Leading Order is determined 
by 



1 + 



Atx 



D 



(33) 



with Uf = 5. Obviously neither the factor 772 nor the redefined operator Oll depend on the 
energy scale fi, the formal dependence within the definition of the operator is compensated by 
the Bag-parameter as given in Eq. 



3.3 Generalization on r]2K 

We want to generalize the considerations leading to an expression of 172 for the i?-system now 
on the i^'-system. Actually this reduces to the question of how to incorporate the effect of the 
thresholds on the renormalization group evolution. In LO this can be seen easily in Eq. (P), 
and a generalization is straightforward. For 772^- we obtain the following result up to NLO 



'n2K 



g(«)/{2/3r) 



{0)x 



9'°V(2/3r) 



W) 



(0)n 



asirrib) 



(Z4 - Z,) + 



w 



D 
'S 



(34) 



So we are now able to derive an expression for the effective Hamiltonian for AS* = 2-transitions 
including QCD corrections up to 0{as)- 



H, 



1 Gl 



eff 



4 7r2 



w 



iv^v*)\ Six,) + iVt,v;fv2 s 

Here, the operator Oll is defined by equation (|3T|) with Uf = 3. 



(35) 



3.4 Gluino mediated corrections 

As should be noted, there are corrections of the various box-diagrams mediated by gluinos as 
well. Basically these corrections exhibit the same topological structure as the gluon-corrections, 
but to some extend they mix the sectors constituting the Flavour Changing Neutral Currents, 
resulting for example in boxes containing W^-bosons, charginos, quarks and squarks (this is 
achieved by replacing the gluon of topology g of Fig. ^). In addition to the non-vanishing mass 
off the gluino this definitely complicates the computational situation considerably. 
Nevertheless, clearly shows, that the vertex-corrections due to the gluinos are of great 
importance and do not vanish in the limit of a heavy gluino-mass. In other words, the gluino 
does not decouple from an effective theory containing only the lighter degrees of freedom. In 
our opinion this indicates that more care has to be taken in dealing with the gluino corrections 
in a consistent fashion. 
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4 Results 



We want to compare now some very well measured experimental quantities related to the 
effective Hamiltonians we investigated. 

BB-mixing can be described by the mass difference AM of the two mass eigenstates, given by 



AM 



|(50|ifeff(A5 = 2)|5)| 



TUB 



M^iyMYr^^SfBrriBBB 



(36) 



where we have employed equation ([l3|) . 

Similar expressions hold for the i^-system. For example, the parameter ex for indirect CP- 
violation in the decay ^ tttt is very well approximated by 



Gp F'^Bkitlk 



M'^Im[V,,V*] 



67r2 y^AMK 

{Re[V,,V:^] [r],S{x,) - r]s{x,, xw)] - Re[Vt,V,:] v^k S} , 



(37) 



under the use of the relation Jm[Krf\4*]* = Im[VtdV^*] [^. 
4.1 Inputs 

Before we present our results, we want to list our input parameters. Note, that we use the 
Wolfenstein-parametrization with parameters A,X,p,ri. We start with the SM-parameters 
and experimental data. They can be found in |5l[ . 



as{Mw) 


= 0.12 , asinib) = 


0.22 , as{mc) = 0.34, 


A 


= 0.80 ±0.08 , A = 


= 0.22, 


mt{Mw) 


= 167 ±10 GeV, 




^BbPb 


= 0.2 ±0.04 GeV 


^Bk/k = 0.135 ± 0.015 GeV , 




= 2.26 ±0.02- 10" 


-3 


AMb 


= 3.1 ±0.1 ■ 10"^^ 


GeV , AMk = 3.5 • 10"^^ GeV 



(38) 

For the MSSM we made, as stated before, some assumptions concerning the parameters. First 
of all, we want to neglect all flavour changing entries in the squark mixing matrix. This leaves 
us with the following form of the two matrices Z entering the final expressions 



( 1 














\ 





1 




















COS0 








sin 











1 




















1 





I 





— sin 








COS0 / 



cos X sm X 
■ sin X cos X 



(39) 
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We dealt with the mixing angles and x free parameters in the ranges [— 7r/4, 7r/4]. We 
further considered for the masses of the squarks and the charginos the ranges 



150 GeV < mi^ < 300 GeV , m^^ < m^-^ < 600 GeV , m^-^ < < 900 GeV , 

100 GeV < m^2, m^, < 400GeV . (40) 

All references on a "reduced" SUSY denote another assumption, namely only one active squark 
and chargino. In other words, in this case we take = x = and xi heavy enough to give no 
more sizeable effects, e.g. decoupled. If not stated otherwise we took the central values of the 
SM-parameters. 



4.2 Results for the 5— system and ck 

It should be noted here, that for the i^^-system there are so far only lower bounds on the size of 
the mass-splitting. They are still somewhat lower than what can be expected from a theoretical 
point of view within the SM. Because the models we investigated here interfere additively with 
the SM we did not explicitly show any figures for the iJ^-system. Of course our results for 772 
and the mass-splitting hold as well, as long as one replaces the flavour d with s when needed. 
Obviously, as can be deduced from Fig.0, the large range for the hadronic parameters of the 
-B-system allows for a large range of the CKM-elements |\4dV^tb| responsible for the mixing 
even in the SM. Furthermore, one can read of the effect of NLO-corrections in comparison 
to the LO result to be of the order of roughly 10%. Fig. § clearly shows the possibly large 
influence of the charged Higgs-sector on the mixing within the 5-system. Obviously this 
sets some constraints on Higgs-parameters within the THDM or even the MSSM, because the 
super symmetric particles interfere additively as well. This is displayed in Fig. where we have 
used our "reduced" MSSM-parameter space. This choice provides a good feeling for the effect 
of the MSSM, by the way, and we can read of possibly large effects on the mixing induced 
by the MSSM. Stated the other way around, the MSSM clearly has the potential to account 
for sizeable effects on the CKM-elements, pushing them to lower values. Note, that similar 
statements hold for the iC-system as well, because the contribution of the c-quarks is not the 
dominant one. 

In the other plots, the scatterplots, we have displayed the effect of NLO-corrections to rj2 as 
compared to the LO-result and found an effect of roughly 10%, the same as in the SM. 
Furthermore, we compared the effect of taking into account the full MSSM using the GIM- 
mechanism for the first two generations within the ranges given above. Clearly, the MSSM has 
the effect of altering the result obtained for the SM. However, as the plots show, the effect of 
a not decoupled charged Higgs is sizeable, and in this sense models like the THDM are worth 
considering. 
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4.3 Effects on the unitarity triangle 

As a matter of fact, the above considerations do influence the investigations concerning the 
unitarity triangle. As shown by some examples (see Fig.|l5|) it will be shifted by adding a 
THDM or a MSSM. 



5 Conclusions 

The investigations above yield a clear picture of the influence of the MSSM on the mixing 
phenomena within the B and the i^'-system. 

Actually it doesn't seem as if NLO-corrections for the MSSM differ considerably from the 
result obtained for the SM. The factor ri2 displayed in (p!OD for our choice of the parameter 
space seems to favour an effect of roughly —10%. However, the NLO corrections do of course 
reduce the uncertainty involved with possibly large scale-differences, but it should be stressed 
that in case these are too large, one has to perform a step-by-step procedure when integrating 
out the heavy degrees of freedom. 



Obviously - as can be seen in Figs.|TT], |I2|, |1^, |TJ - the MSSM might have some considerable 



influence on the observables under consideration and it has clearly the potential to shift the 
point (p, rj) of the unitarity triangle in the Wolfenstein parametrization. However, an investi- 
gation of the effect of gluinos and neutralinos on the NLO-level is still missing and this should 
be done to have a flavour of the complete picture as predicted by the MSSM. Furthermore, one 
should investigate the influence of gluino-mediated QCD corrections as accomplished for the 
decay 6 57 by [pl|]. 

Nevertheless we have to remark, that the large range of the Bag-parameters shadows the effect 
of physics beyond the Standard Model. 
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A Inami— Lim Functions 

We consider only the very region of low tan/3 1. No scalar operator emerges for the effective 
Hamiltonian. 



A.l Leading Order Inami— Lim functions 

bwwKXw) jf. prr; 7^ , 
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A.2.3 MSSM, chargino sector 
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Figure 1: Diagrams contributing to BB-mixing in LO within the SM. The would- 
be Goldstone bosons $ emerging when using the Feynman-t'Hooft gauge are displayed 
explicitly, crossed diagrams are missing. 




Figure 2: Diagrams contributing to SB-mixing in LO 
within the THDM. The would-be Goldstone bosons $ are 
subsumed in the labels W, crossed diagrams are missing. 
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Figure 3: Diagrams contributing to BS-mixing in LO witiiin the MSSM. Only the 
boxes involving supersymmetric particles are displayed. Note that the neutralinos and 
gluinos arc Majorana-particlcs. This fact produces graphs like the last one displayed 
where we have shown explicitly the effect of the Majorana's character by adding the dots 
and the direction of the flow with arrows. This feature results in scalar operators. 
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Figure 4: Diagrams contributing to the factor r/2 for the QCD-corrections in LO. 
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Table 1: The Fcynman rules employed throughout this article. Note, that only the 
part proportional to the strong coupling constant enters the four-squark vertex. 
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Figure 6: Diagrams responsible for explicit QCD-corrections contained in 772 at NLO 
within the chargino-squark sector of the MSSM. 
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Figure 7: The product of the CKM-elements 
l^d^bl versus ^/BbTb allowed region of rrit 

within the SM. 




Figure 8: The influence of a pure THD M for d if- 

fcrcnt values of Mh and tan/3. We chose Bsfg = 
0.18 GeV, \VtdV;f,\ = 0.007 and mt{Mw) = 167 GeV. 
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Figure 9: The ratio of the mass-splittings AM 

obtained for our "reduced" MSSM and the SM. 
We assumed a Higgs-mass of Mh = 500 GeV and 
tan/3 = 1.5. 
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Figure 10: The ratio ri2Lo/V2NLO in dependence 
on rrif^, "^x^^ X- In cases we practically decou- 
pled the Higgs boson by setting Mh = 1000 GeV and 
tan/3 = 5. 
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Figure 11: The ratio of the mass-spHttings AM 
within the MSSM and the SM in dependence on rrif^, 
m^2, X- In S'll cases we practically decoupled the 
Higgs boson by setting Mh — 1000 GeV and tan/3 = 
5. 
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Figure 12: The ratio of the mass-spHttmgs AM 
within the MSSM and the SM in dependence on m^^, 
m^^ , X- III £^11 cases we included the Higgs boson by 
setting Mh = 100 GeV and tan/3 = 1.5. 
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Figure 13: The ratio of ex within the MSSM 
and the SM in dependence on rrif^, ^x2j X- In 
all cases we decoupled the Higgs boson by setting 
Mh = 1000 GeV and tan/3 = 5. 
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Figure 14: The ratio of ex within the MSSM 
and the SM in dependence on m^^, m^^, x- In 
all cases we included the Higgs boson by setting 
Mh = 100 GeV and tan/3 = 1.5. 
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Figure 15: The influence of adding THDM and 

MSSM to the SM on the position of (p, 77). We varied 
the top-mass within the bounds indicated above as 
well as the hadronic parameters and chose Mh = 
200 GeV, tan/? = 1.5 for the THDM investigations 
and Mh = 500 GeV, tan/3 = 5, = 200 GeV and 
m^, = 120 GeV for the reduced MSSM. 
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